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Abstract 

In our calculation of neutron star crust heating we include several key new model features. 
In earlier work electron capture (EC) only allowed neutron emission from the daughter ground- 
state; here we calculate, in a deformed QRPA model, EC decay rates to all states in the daughter 
that are allowed by Gamow- Teller selection rules and energetics. The subsequent branching ratios 
between the In,. . . ,xn channels and the competing 7-decay are calculated in a Hauser-Feshbach 
model. Since EC accesses excited states, many more neutrons are emitted in our calculation than 
in previous work, leading to accelerated reaction flows. In our multi-component plasma model a 
single (EC,xn) reaction step can produce several neutron-deficient nuclei, each of which can further 
decay by (EC,xn). Hence, the neutron emission occurs more continuously with increasing depth 
as compared to that in a one-component plasma model. 
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Nuclear reactions in the crust of a neutron star (NS) form an important bridge between 
the extensively studied nuclear burning on the surface and the less well understood neutrino 
cooling processes in the core, where exotic states of matter can exist. Observation of ac- 
creting neutron stars in binary systems is proving to be a promising new tool to probe the 
thermal evolution of the inner core and its neutrino cooling rate. In particular, the ignition 
of "X-ray superbursts" (XRSB) due to ^^C fusion has been shown to be sensitive to the 
thermal profile of the neutron star crust Similarly, neutron stars exhibiting intermittent 
periods of rapid accretion ( "transients" ) have a quiescent luminosity which is thought to be 
powered by the heat released in the crust 

Surface burning on an accreting NS by the "rp-process" (rapid proton capture) in a Type- 
1 X-ray Burst (XRB) produces heavy elements near the proton drip-line. As these 

ashes are buried by subsequent accretion successive electron capture (EC) reactions occur as 
the electron chemical potential E-p increases rapidly with depth. In pioneering work a decade 
ago, Haensel and Zdunik outlined these EC reactions in a one-component plasma (OCP) and 
found that the heat released in EC was not substantial - most of the energy release occurred 
much deeper in the crust through pycnonuclear fusion reactions . More recently in a 
multi-component plasma (MCP) simulation, Gupta et.al. found that EC to excited states 
and subsequent 7-deexcitation increased the heating at shallow depths j?]. In this letter 
we report that EC captures to highly excited daughter states also lead to neutron emission 
processes near the neutron drip-line that are vastly accelerated when the level schemes of 
deformed neutron-rich nuclei are incorporated as opposed to when emissions are restricted 
to ground states (g.s.). These reactions fundamentally change the direction and speed of 
the reaction pathway and alter the crustal heating profile. 

For Ep > 27 MeV there will be a neutron chemical potential /in ^ and consequently 
a free neutron abundance. The depth Ep ^ 27 MeV is called the "neutron drip" (ND) 
point. Close to ND {Ep ~ 25 MeV), EC can occur into excited states with Ecxc > •S'xn 
where is the separation energy for x neutrons. These reactions are henceforth referred 
to as (EC,xn) reactions. Once (EC,xn) is allowed at Ep ~ 25 MeV, any EC is always 
accompanied by some neutron emission since the x ^ 3 channels are always open and the 
X = channel ("pure EC") is highly suppressed. Post-ND the products of (EC,xn) can be 
highly EC- unstable ( "superthreshold" ) at [Z — 1) since the neutron number has also been 
substantially lowered. Although some of these nuclei will recapture neutrons, yet the favored 
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after (n,7) will not be as high as the pure EC pathway would have permitted, due to the 
powerful non-thermal driving force removing neutrons. This force is [Ep — -Ethresh.g.s. — •S'xn) 
and since it is much greater than kT the superthreshold (EC,xn) rates are also independent 
of temperature. The temperature can only drive neutron emission for S^^ ^1.4 MeV by 
(7,n). Consequently, the density-driven (EC,xn) rates, which are more than a factor of 
10^ larger than the (7,n) rates, expand the path into a band 1.0 MeV < 5*11 < 3.5 MeV. 
Nuclei in this band are generally EC-unstable at Ep ^ 25 MeV. Hence, each of the (EC,xn) 
products (resultant network nodes for different x) are subjected to reaction-driving energies 
{Ep — Eth_Tcsh,g.s. — 'S'xn) that can be higher than in the previous (EC,xn) reaction step, leading 
to a dense web of subsequent (EC,xn) branchings or a "Superthreshold (EC,xn) Cascade" 
(henceforth called the SEC-process) rapidly lowering both (Z) and (A). For clarity, we 
denote the absolute value of the g.s. nuclear mass differences as -Ethrcsh,g.s. while an actually 
encountered threshold is always written as -Ethrcsh and includes the energy of the first state 
to which EC is allowed, which may not be the g.s. 

Our network includes EC (/?"'") and /5~ rates from the proton drip-line to the neutron drip- 
line using transition matrix elements from the QRPA model described in l8|. Calculated 
values of deformations from [9| were consistently incorporated for Z > 8 nuclei , below 
which we have assumed a typical quadrupole deformation 62 = 0.55. From any excitation 
energy E^xc in the EC daughter {Z — 1,A) the neutron decay and competing 7-ray emission 
branching ratios are calculated with the nuclear reaction code GNASH 10|]. This code 
employs the statistical Hauser-Feshbach (HF) method augmented with a Direct Semi-Direct 
(DSD) component for those nuclei in which the direct reactions will dominate. The 
maximum excitation energy considered is 30 MeV, and up to 20 emitted neutrons are allowed. 
We account for neutron kinetic energies which reduce the available energy for subsequent 
neutron emissions. For Z < 10 we do not use HF or DSD for the neutron branches: rather, 
if EC occurs into a state above S'xn but below S'(x+i)n, then the x-neutron branch is given 
full strength and all others are set to zero. 



We use radiative neutron capture (n,7) rates from [12|]. The reverse (7,n) rate for a 
typical pre-ND crust nucleus ^^^^Kr (ignoring partition functions) is A = {Nx{av)Y^) x 
10^° X Tg^^ exp{-ll.Q05S,,/Tg) = 2.74 x 10~S~^ for pre-ND neutron abundance = 10-^ 
a forward (n,7) cross-section Nx{av) = 3 x 10^ cm'^g~^s~^ and Sn = 1.348 MeV at crust 
temperature 0.5 GK. Moving towards a higher = 2.079 MeV in ^''^Kr with NAicrv) = 
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2 X 10^cm^g~^s~^, we have A = 7.82 x 10~^^s~^ which is aheady too slow to compete 
with the accretion timescale At ~ lO^s during which the resuhing AEp can make (EC,xn) 
reactions competitive. Typically, only nuclei with ^ 1.4 MeV have (7,n) competitive 
with accretion at T ~ 0.5 GK. The products of these (7,n) are susceptible to EC which 
produce {Z — 1) nuclei, and the process continues to lower A by further (7,n). In models 
that rely on the stellar photon bath to produce neutrons by (7,n) reactions, this will be a 
valid pathway pre-ND (but constrained by ^ 1.4 MeV). In the vicinity of ND [Ep ~ 25 
MeV), however, the (EC,xn) rates are more than a factor of 10^ larger than the (7,n) 
rates. Also, from a heating efficiency perspective, the slow pre-ND (7,n) reactions have a 
negligible effect - beyond Ep ~ 25 MeV the SEC-process removes any residual compositional 
memory. Furthermore, at Ep > 20 MeV the number density of thermal photons is suppressed 
because the propagating photon is dressed in the medium (the dressed particle is called the 
plasmon) and acquires a finite mass htOp = {AnNee'^ /me){l + (/i/(mec))^(37r^Ai'e)^/^) ^ 1.3 
MeV at pYe ~ 10^"*^ g/cc increasing to about 4.1 MeV at 10^^ g/cc. Here Yg, TVg, me are 
electron fraction, number density and rest mass, respectively. Under typical conditions in 
the crust (10^^ g/cc, 0.5 GK) the number density of photons is suppressed by the factor 
exp{—hLjp/kT) ~ 10""*^^. Since hoUp > ~ 1.4 MeV, this suppression has a strong effect 
at energies important to the pre-ND (7,n) reactions. Furthermore, the plasmon suppression 
increases with increasing density: exp{—hujp/kT) ~ 10~^^ at 10^^ g/cc and 0.5 GK. However, 
the higher the electron density, the more efficient the (EC,xn) processes become at removing 
neutrons from increasingly higher excitation energies. 

Furthermore, the high Nuclear Level Densities (NLDs) at Ecxc ~ 'S'xn (in-d the high su- 
perthreshold phase-space factors involved ensure that the SEC-process is a robust mechanism 
with negligible sensitivity to uncertainties in the underlying nuclear model. The (EC,xn) 
timescales will always be shorter than the accretion and the (7,n) timescales by several or- 
ders of magnitude in a global nuclear model used to calculate the EC daughter excitation 
energies. The speed of the resulting reaction fiow will be unaffected by the choice of global 
nuclear model, though the "transit points" at a given Ep may differ slightly. We highlight 
two aspects of the QRPA model we have used which have significant global impact: (1) the 
single-particle model serving as input to the QRPA is a complete model space (no levels will 
be systematically missing) and, (2) we treat deformation realistically. Though global Shell 
Model calculations of weak interaction strengths for the (heavy and/or neutron-rich) nuclei 



4 



25 26 27 28 29 30 31 32 33 34 35 
Electron Chemical Potential £f (MeV) 

FIG. 1: (color) Abundance distributions Y(N) and Y{Z) shown for abs(log;^Q(y)) < 5 at AE-p = 0.1 
intervals in the E-p range 25-35 MeV for an MCP calculation starting with pure ^''^Pd. Beyond 
Ep ~ 26 the MCP distributions of Y^N) and Y{Z) rapidly broaden due to EC into highly excited 
states followed by neutron emissions. The OCP trajectories of [l3| starting with the same ICC and 
evolving through EC and g.s. neutron emissions lie outside the 2a limits of the MCP Y{N) and 
Y{Z) distributions. 

of interest are simply not available, the same considerations of high NLD at high i^exc apply 
to nuclear structure computed in the Shell Model as they do to the QRPA calculations. 
The combined effect of g.s. to g^. EC and g.s. neutron emissions in the crust was referred 



to as the "inverse r-process" by [ij]. However, the possibility of "superthreshold" (EC,xn), 
with very large driving forces that lead to further cascades within a narrow Ai^p slice of the 
crust was not explored. In fact, since earlier models had neutron emission only from g.s. 
they missed important nuclear structure effects in deformed neutron-rich nuclei. The level 
spacings in these exotic nuclei can influence flows by changing (EC,xn) rates and branchings 
with only a small change in Ep. Moreover, (1) the SEC-process is quite distinct from the 
r-process in the Waiting Point Approximation where at a certain temperature (n, 7) — (7, n) 
equilibrium along constant Z occurs since the (EC,xn) "reverse rate" is more than a factor 
of 10^ larger than the (7,n) rate (even without plasmon suppression); (2) the (EC,xn) rate 
increases with proton-richness, unlike the (7,n) rate in the r-process; and (3) the non-thermal 
reaction energy source Ep is not even approximately constant, as it increases steadily with 
depth in the NS crust. Therefore, there are few similarities with the r-process, except that 
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that the SEC-process runs in the opposite direction in the ZA^-plane. 

Since the SEC-process covers wide ranges of the nuclear chart with a dense web of compet- 
ing (EC,xn) channels, it becomes difficult to identify individual nuclei that play a dominant 
role at a given Ep once the r.m.s. deviation AA^ increases to ^ 5 for the abundance distri- 
bution ^(A^). Rather, we plot the means (A^) and (Z) and the r.m.s. deviations AZ, AN 
(and also 2AN and 2AZ) of the abundance distributions Y{N) and YiZ) with changing 
E-p in fig{T]to compare our MCP evolution with the OCP evolution of 13|] starting with the 
same initial crust composition (ICC) of pure ^°^Pd. We start the evolution at Ep 0.5 MeV 
(p ~ lO^g/cc) and terminate our simulation when fusion reactions begin to cycle material 
out of the Z ~ 10 region at Ep ~ 38 MeV (p ~ lO^^-Sg/^^)^ 

Prior to neutron emission the "pure EC" path through even-even nuclei is the same as in 
the OCP calculation of i06pd^i°6Ru-*i°6Mo^i°6Zr^i°6gr^i06j^r^i06gg_ However, 



i06Sg_^i06Ag at Ep = 25.86 MeV is to Ecxc = 2.45 MeV > S3,, = 0.99 MeV and neutron re- 
moval can produce nuclei up to ^°^As. Now ^thresh = 23.46, 24.44, 24.12 MeV for i03.i04,i05^g 
respectively allow "superthreshold" (EC,xn) branching with high values of x. Thus, (A^) 
begins a rapid descent and the MCP trajectory does not produce A^ = 74 at ^''^Ge accessed 
at Ep ~ 27 MeV in the OCP (fig.[I]). Well-deformed ^"^As (eg = 0.225 from Q) has neu- 
tron energy levels spaced by a few 100 keV, each participating in (EC,xn) over a very small 
change AEp in Ep. This effect combined with a low S,, results in multiple neutron emission 
towards (spherical) subshell closures at A^ = 64, 56. Hence, AA^ suddenly expands to ~ 5 
beyond the subshell A^ = 70 and the slope of (A^) abruptly changes as neutrons are rapidly 
emitted between A^ = 66 and 56. We follow the increasingly dense web of (EC,xn) branch- 
ings from ^°^As(EC,ln)^'^^Ge onwards. The product undergoes up to (EC,3n) at Ep ~ 28 
MeV. A dominant branch at 53% is to ^°^Ga whose (EC,xn) products are ^°°~^^Zn, each 
with branchings between 30 and 40%. The flow is now no longer concentrated along a single 
x-neutron emission channel. At -Ethrcsh = 28.6 MeV ^^Zn(EC,xn) occurs with x < 3. The 
product ^^Cu (A^ = 66) has -Ethrosh = 28.8 MeV and releases up to 6 neutrons with sub- 
stantial branching ratios. Due to this rapidly expanding web of (EC,xn) with increasing x, 
neutrons are emitted continuously as Ep rises, and consequently at Ep ~ 29 MeV (fig. [I]) 
fiow from A^ = 66 to A^ = 60 is already at the factor of 3 level by abundance. In contrast, 
the OCP essentially does not evolve from Ep = 27.08 up to 31.22 MeV, i.e. there is no ac- 
tivity between production of ^°^Ge and its destruction by ^'^^Ge-|-4e~-^^^Ni+14n+4z/e which 
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abruptly releases 14 neutrons at a single Ep = 31.22 MeV. This scenario is highly unlikely 
if (1) neutron emissions occur from excited states - instead of the typical (2EC,6n) steps 
of the OCP evolution at widely separated Ep, (EC,xn) with much higher x can participate 
in the MCP at several closely spaced Ep, and (2) there is more than one network node 
so that 100% of the flow is not controlled by a single x ~ 3 neutron emission mode like 
a typical OCP (2EC,6n) on *'^Cr at Ep = 32.76 MeV. If we remove only 3 neutrons after 
the threshold capture on '^''Cr, the product ^^V is still superthreshold to (EC,xn) by almost 
4 MeV. The MCP branchings ^ 10% on the ^''V intermediate nucleus are 12.6%(EC,4n), 
38.4%(EC,5n), 15.8%(EC,6n), 19.2%(EC,7n) showing odd-even and neutron kinetic energy 
subtraction effects. The resulting large spreads AA^, AZ in the MCP evolution allow much 
smoother movement through several MCP network nodes. This occurs because any one of 
the network nodes [Z — 1, N + 1), [Z — 1, N), [Z — 1, N — 1), ... resulting from different x in 
a single (EC,xn) step can itself move on to lower and Z through (EC,xn) either at the 
prevailing Ep or with a small increase AEp ~ 0.3 MeV. 

The nuclear structure probed by the high energy {Ep — -Ethresh) of "superthreshold 
(EC,xn)" determines the characteristics of the reaction flow in the following ways: 

{l)The slopes of {N),{Z) with Ep, or the speed of the reaction flow, comes from 
a removal of the (spherical) subshell degeneracy due to deformation. This effect is 
strongest midway between spherical (suh)shells. The OCP curve passes through = 
64,60,56,52,48,44,40,36,30 with large gaps AEp ~ 1 MeV(and a huge gap of 4.14 MeV 
between the productions of ^°^Ge and ^^Ni). In contrast, the MCP neutron emission is much 
more continuous with changing Ep. At higher Ep in the MCP evolution, fig. [1] clearly shows 
that within the (spherical) neutron ggj2 subshell (A^ between 40 and 50) there are several 
transit points in A^ that become accessible at very closely spaced Ep. These correspond to 
the well- separated deformed neutron levels emanating from the same (spherical) subshell. 
A level diagram of mid-proton-shell *°Cr (deformed with e2 = 0.233) not only shows this 
effect, but also that upon entry into the A^ = 56 subshell closure, A^ = 54, 52 should play a 
prominent role and A^ = 52, 50 correspond to levels at almost identical energy. Therefore, 
the MCP A^-evolution flattens at A^ = 56 and 40, but not at A^ = 50, through which it is 
rather steep. Throughout the evolution the OCP slope of N{Ep) (with notable exceptions 
at A^ = 64,60 corresponding to Z = 28,26) is closer to the MCP slopes of (A^)(±2AA^) 
between major neutron (sub)shells (see the MCP slopes at A^ = 56,40 in fig.[I]) - thus the 
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OCP g.s. neutron removals do not reflect the rich structure within neutron (sub)shells for 
well-deformed neutron-rich nuclei. 

(2) AA^ and increase to reflect the occupancies of accessible (suh)shells at a given Ep. 
An example is the increase of AA^ to (40 — 28)/2 = 6 closely matching the network value of 
6.44 at Ep = 33.5 MeV. Levels beneath subshell closures at A^ = 38, 32 are all accessible at 
this Ep but A^ between 28 and 20 is only accessed at A^ = 26 by highly abundant nuclei. 

(3) Unlike the OCP evolution, branchings populate odd-A chains also - these are rapidly 
depleted, driving abundance to lower Z faster. 

(4) The free neutron abundance increases smoothly with Ep due to the large number 
of available (EC,xn) sources at each depth in marked contrast to the discontinuous evolution 
ofX^{Ep) in the OCP (see fl, Q, Q/;. 

The resulting distribution ^(A^) at any Ep is multiple-peaked primarily due to odd-even 
and neutron (sub) shell effects. It is also highly asymmetric with long neutron-deflcient tails 
that are susceptible to further (EC,xn). This leads to a broadening ofY{Z) to the extent 
allowed by the accessible proton (sub)shells. 

Beyond Ep ~ 29 MeV the MCP rapidly diverges from the OCP since roughly [(2AZ') x 
(2AA^)] ^ 5 X 10 = 50 network nodes are simultaneously processing the abundances unlike a 
single node in an OCP scenario. This dense Markov process of (EC,xn) transitions (1) quickly 
removes compositional memory of the XRB ashes through several A^-nodes at each Z (and 
vice- versa), which is not possible in an OCP evolution, and (2) changes the evolutionary 
description from one of individual nuclei at speciflc Ep (which is possible for the OCP) into 
that of r.m.s. deviations and higher moments which ultimately determine the shapes of the 
y (A^) and Y{Z) distributions. These shapes are clearly discernible in flg.[T] since the color 
changes denote abundance variations by a factor of only lO"'^ ~ 3. The large spreads AA^ 
and AZ in these distributions are important inputs to models of transport processes in the 
crust such as those determining thermal conductivity. 

The efficiency of the SEC-process in heating the crust is shown in fig. El We have the re- 
markable result that MCP crustal heating profiles from two very different ICC from possible 
Superburst progenitors change slope drastically near Ep ~ 35 MeV and converge to within 
~ 0.1 MeV/u of each other. This is of the order of the difference in BE/A of the ICC. Note 
that the MCP heating diverges from the OCP at a depth Ep ^ 30 MeV where the MCP 
slopes of (A^)(±2AA^) change dramatically in fig.[Il 
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FIG. 2: (color) Crustal heating in our MCP model for ICC (XRB ashes from Q]) from NS accreting 
5 % solar (x) and solar metallicities {X) at the accretion rate M = O.lmEdd characteristic of Su- 
perburst progenitors (accretion at the Eddington limit is rriEdd ~ 8.85- 10^gcm~^s~'^). Also plotted 
are OCP heating from [2] for the ICC comprised of ^'^^Pd. ICC from multizone 

XRB model xM: Peak A ~ 100 — 110 (very similar to the 1-zone rp-process products of [j] and 
also to an ICC of pure ^^^^Pd in the OCP model of Q])- ICC from model XM: A 

bimodal abundance distribution with significant pre-Fe-peak abundance, and also in yl ~ 60 — 80. 

The bulk of deep crustal heating by the SEC-process occurs in a narrow region post-ND, 
rather than near the crust-core interface if pycnonuclear fusion dominates. Being further 
from the core this new crustal heating mechanism can result in thermal profiles conducive 
to the ^^C ignition of X-ray Superbursts. Furthermore, unless the ICC is comprised only of 
pre-Fe-peak nuclei concentrated in A = 36 — 44 (fig. 7 of j^j explains why shallow heating 
should differ in this case; however the total crust heating will still differ by ^ 0.2 MeV/u from 
MCP evolution of other ICC), the composition as a function of depth is not dependent on 
the ICC. Thus for a large class of accreting NS, the composition-dependent conductivities 
and crust neutrino cooling will look very similar, and so will the crustal thermal profiles 
which depend on these quantities. 
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